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Abstract 

Using the full potential linear muffin-tin orbitals (FP-LMTO) method we examine the pressure- 
dependence of superconductivity in the two metallic phases of Boron: bet and fee. Linear response 
calculations are carried out to examine the phonon frequencies and electron-phonon coupling for 
various lattice parameters, and superconducting transition temperatures are obtained from the 
isotropic Eliashberg equation. The fee phase is found to be stable only at very high pressure 
(volume per atom < 21.3 bohr^), estimated to be in excess of 360 GPa. The bet phase (volume 
per atom > 21.3 bohr^) is stable at lower pressures in the range of 210-360 GPa. In both bet 
and fee phases the superconducting transition temperature Tc is found to decrease with increasing 
pressure, due to stiffening of phonons with an accompanying decrease in electron-phonon coupling. 
This is in contrast to a recent report, where Tc is found to increase with pressure. Even more 
drastic is the difference between the measured Tc, in the range 4-11 K, and the calculated values 
for both bet and fee phases, in the range 60-100 K. The calculation reveals that the transition 
from the fee to bet phase, as a result of increasing volume or decreasing pressure, is caused by 
the softening of the X-point transverse phonons. This phonon softening also causes large electron- 
phonon coupling for high volumes in the fee phase, resulting in coupling constants in excess of 2.5 
and Tc nearing 100 K. 

Although it is possible that the method used somewhat overestimates the electron-phonon cou- 
pling, its success in studying several other systems, including MgB2, clearly suggests that the 
experimental work should be reinvestigated. We discuss possible causes as to why the experiment 
might have revealed Tc's much lower than what is suggested by the present study. The main 
assertion of this paper is that the possibility of high Tc, in excess of 50 K, in high pressure pure 
metallic phases of boron cannot be ruled out, thus pointing to (substantiating) the need for fur- 
ther experimental investigations of the superconducting properties of high pressure pure phases of 
boron. 

PACS numbers: 74.70.-b, 74.62.Fj, 74.25.Kc, 71.20.-b 
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I. INTRODUCTION 



Although metalhzation of boron under pressure was predicted on the basis of electronic 
structure calculations some time agoi , experimental verification of this result has so far 
been lacking. On the basis of the calculations^ nonmetallic icosahedral boron(B) is expected 
to undergo a structural transition, first to a body-centered tetragonal (bet) phase at ~210 
GPa and then to a face-centered cubic (fee) phase at ~360 GPa. Recently high pressure 
experiments by Eremets et al^ have found boron to be not only a metal, but also a super- 
conductor at high pressure, with superconducting transition temperature Tc increasing with 
increasing pressure. Superconductivity appears at around 160 GPa, and increases from 
6 K at 175 GPa to 11.2 K at 250 GPa. The appearance of superconductivity in B under 
pressure is not surprising, because the existence of different phases is usually an indication 
of strong electron-phonon coupling, which is (most often) also responsible for superconduc- 
tivity. In this work we study the superconductivity in bet and fee phases of B using ab initio 
theoretical calculations. We use the linear response scheme implemented by Savrasov^i^ to 
study the phonon properties and electron-phonon coupling. The electronic structure is cal- 
culated by using the full potential linear muffin-tin orbital (FP-LMTO) method^. Recently 
such calculations have been performed for MgB2^ and several other systems^i^ including 
some boron containing compounds^'^°, and have been found highly successful in capturing 
the essence of electron-phonon superconductivity in these compounds. 

A theoretical study of superconductivity in the fee phase of boron was reported by Pa- 
paconstantopoulos and Mehlii shortly after Eremets et al? reported their experimental 
results. However, the lattice parameters for the fee phase used by the authors correspond 
to volumes per atom for which a lower energy structure is bet, as had been revealed by 
an earlier calculation of Mailhoit et al}. Papaconstantopoulos and Mehlii calculate the 
electron-phonon coupling or the Hopfield parameter using the rigid muffin-tin approxima- 
tion, arguing that the fee phase might be metastable for such lattice parameters. This 
speculation is based on their calculation of elastic constants showing Cu — C12 > 0; Cu > 0, 
i.e., the Born criteria — for the stability of the cubic phase are not violated. However, our 
linear response calculations yield imaginary phonon frequencies for lattice parameters in 
excess of 4.4 a.u. in the fee phase. Thus, although the existence of a metastable fee phase 
for lattice parameters 4.6-6.0 a.u. used by Papaconstantopoulos and Mehl^i cannot be ruled 
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out, the chances are that the equihbrium phase at such lattice parameters is not fee. Our 
calculations for equilibrium fee phase (lattice parameter less than 4.4 a.u.) shows a volume 
dependence of the superconducting transition temperature that is the opposite of what is 
observed experimentally, i.e., decreases with increasing pressure. The study by Papacon- 
stantopoulos and Mehlii only involves the fee phase with large lattice parameters, for which 
the structure is unstable or, at best, metastable. The present work studies both the fee 
and bet phases, each within its appropriate range of lattice parameters. Our linear response 
calculations also indicate the onset of a bet phase via preference for a change in the c/a ratio 
as the volume per atom is increased. Detailed discussion of the relevance of our results to the 
experimental work and that of Papaconstantopoulos and Mehiii is provided in appropriate 
sections. 

The organization of this paper is as follows. In section II we discuss the stability of 
the bet and fee phases of boron as a function of volume per atom. In Sections III and IV 
we discuss the calculated electronic bands, phonon spectra, and electron-phonon coupling 
for the fee and bet phases, respectively. In section V we discuss the volume-dependence of 
the superconducting transition temperature, summarize and compare our results with the 
experimental data and discuss possible sources of difference between the two. 

II. STABILITY OF FCC AND BCT PHASES 

The volume vs. energy curves for fee and bet phases of boron, obtained via FP-LMTO 
calculations, are shown in Fig^ A constant energy of 49.0 Ry. has been added to the 
energies per atom for convenience in plotting. These results are very similar to those given by 
LMTO-ASA method, and are in close agreement with the earlier plane-wave pseudopotential 
calculations of Mailhoit et al}. Pseudopotential calculations of total energy of the bet 
structure by Mailhoit et al} as a function of c/a had indicated a global minimum around 
c/a 0.6. Via an elaborate study of the total energies for various monoclinic and tetragonal 
distortions of the fee unit cell as well as the bet total energies for various c/a values these 
authors concluded that the bet minimum occurs around c/a ^ 0.65. We have not carried 
out extensive calculations to locate the exact value of c/a, but our FP-LMTO calculations 
on a crude mesh of c/a for fixed volumes do indicate a value in the range 0.6 — 0.7. Our 
linear response calculations for c/a = 0.65, however, produced complex phonon frequencies 
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FIG. 1: Energy vs. volume (per atom) curves for the foe and bet phases of boron. 

at the symmetry point N and nearby wave vectors. Increasing c/a to 0.675 resulted in real 
frequencies at all symmetry and intermediate points. Thus the linear response results in 
this paper are presented for this value of c/a. It is evident from Fig^that the actual energy 
values for c/a = 0.65 and 0.675 for the same volume in the bet phase are very close, and 
their energy- volume curve intersects the fee energy- volume curve at almost the same point. 
Thus the conclusions regarding the stability of the bet phase against fee are not influenced 
critically by the choice of c/a in the range 0.65 — 0.675. 



III. THE FCC PHASE 
A. Energy bands 

In FigI21 we show the FP-LMTO energy bands in fee B as a function of the lattice 
parameter 'a'. Calculations were carried out with a double-K spd LMTO basis set for 
describing the valence bands. The charge densities and potentials were represented by 
spherical harmonics with / <6 inside the nonoverlapping (touching) muffin-tin spheres and 
by plane waves with energies <495 Ry in the interstitial region. Brillouin zone integrations 
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were carried out with the tetrahedron methodic using a 24,24,24 division of the Brillouin 
zone, corresponding to 413 wave vectors in the irreducible part. The extraordinarily large 
band- widths in FigO are due to the extreme pressures of the systems. The fee phase is 
stable for a=4.4 a.u. and less, and unstable for a=4.6 a.u. and higher values. The smallest 
volume per atom at which the bet phase is stable is around 3.55 A'^, according to the 
volume-energy curves shown above. This volume corresponds to a fee lattice parameter of 
4.576 a.u. One of the two energy bands crossing the Fermi level between the two symmetry 
points X and W becomes very flat at larger volumes. The character of the fiat energy band 
between X and W is primarily 's' and 'p': at W it is purely of 'p' and at X it is purely 
of 's' character. Papaconstantopoulos and MehUi have speculated that this fiat band is 
responsible for superconductivity of the fee phase in boron, pointing out that a similar fiat 
band at the Fermi level appears in MgB2. Fermi surface plot of Papaconstantopoulos and 
Mehiii shows hole pockets near X and electron poclcets near W, as would also be expected 
according to the present calculations. Although we have not verified this explicitly, it is very 
likely that the fiat band between X and W points contributes significantly to the electron- 
phonon coupling in fee boron. As this band becomes flatter for larger lattice parameters, 
the electron-phonon coupling increases. However, this large electron-phonon coupling also 
makes the fee phase unstable at higher volumes via softening of X-point transverse phonons, 
and brings about the fcc^ bet transition, as discussed below. The steep band crossing the 
Fermi level between X and W is of 'p' character. 

B. Phonons and Electron-phonon Coupling 

In Fig|21we show the variation of the calculated phonon spectra with the lattice parame- 
ter for fee B. The phonon frequencies were calculated for 29 wave vectors in the irreducible 
part of the Brillouin zone (BZ), resulting from a 8,8,8 division of the BZ. This division of 
the BZ yields only a small number (3-4) of wave vectors along the symmetry directions. 
The solid lines in FigOl should be taken only as a guide to the eye, rather than actual 
phonon branches with correct branch crossings. With increasing volume, the phonon fre- 
quencies soften throughout the Brillouin zone, as expected. However, the softening of the 
two transverse phonon branches near the X point is most pronounced. The X-point trans- 
verse phonons are related to the 180" out of phase vibration of the two atoms in a bet 
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FIG. 2: FP-LMTO energy bands in fee Fe for four different lattice parameters. The flat energy 
band between the symmetry points X and W, which provides strong electron-phonon coupling, also 
renders the structure unstable at higher volumes. The zero of energy has been set at the Fermi 
level. 

unit cell (c/a= 1.414 for fee lattiee). Softening of these plionons with increasing volume 
indicates growing instability of the fee structure with respect to the c/a ratio, and acts as 
a precursor to the fcc^bct phase transition. This phonon-softening has two important ef- 
fects: it increases the phonon density of states at low frequencies beyond the usual parabolic 
DOS given by the Debye (continuum) model and increases the electron-phonon coupling in 
the low frequency region (partly due to increased phonon DOS and partly due to increased 
electron-phonon matrix element). These effects are clearly seen in FigjH where the calcu- 
lated Eliashberg spectral function (upper panel (a)) and the phonon density of states (lower 
panel (b)) are displayed for three different lattice parameters. The sharp peaks in both the 
density of states and the Ehashberg spectral function at the high end of the spectrum are 
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FIG. 3: Calculated phonon spectra for foe B for the lattice parameters (a) a = 4.0 a.u., (b) a = 
4.4 a.u. and (c) a = 4.6 a.u. The dots represent the calculated phonon frequencies, with solid lines 
providing only a guide to the eye. The softening of the transverse phonons at the symmetry point 
X occurs as the system expands to a volume above 3.4 A'^/atom, in accordance with the earlier 
results of Mailhoit et al. (RefO). 

due to the lack of dispersion in the longitudinal phonon band close to the X point in the 
F-X direction. The lack of dispersion is enhanced as the atoms move further apart, giving 
a higher density of states and consequently higher electron-phonon coupling. 

In short, the transition from fee to bet phase with increasing volume (lowering of pres- 
sure) is driven by softening of transverse X-point phonons, driving the system to a c/a 
ratio different from (in this case, lower than) the fee value of 1.414. The system in this 
pressure/ volume region shows strong electron-phonon coupling, suggesting the possibility of 
superconducting phase with a relatively high transition temperature Tc. 
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FIG. 4: Eliashberg spectral function (a) and phonon density of states (b) for fee B, for three 
different lattice parameters. 

IV. THE BCT PHASE 
A. Energy bands 

In Fig0we show the FP-LMTO bands structure of bet B for a=4.35 a.u. for two different 
c/a ratios: 0.675 and 0.65. Another band structure for c/a=0.675, but a= 4.0 a.u. is also 
shown for comparison. Calculations were carried out with a double-/? spd LMTO basis set 
for describing the valence bands. The charge densities and potentials were represented by 
spherical harmonics with / <6 inside the nonoverlapping (touching) muffin-tin spheres and 
by plane waves with energies <400 Ry in the interstitial region. Brillouin zone integrations 
were carried out with the tetrahedron method^? using a 30,30,30 division of the Brillouin 
zone, corresponding to 1992 wave vectors in the irreducible part. 

As in the fee case, the band widths are abnormally large due to extreme high pressures 
of the systems. The change in the band structure is understandable as a result of decreasing 
volume from (a) to (b). In panel (c) of FigElwe show the bands for a c/a ratio of 0.65, for 
which one of the transverse N-point phonons has imaginary frequency. There is no discern- 
able feature in the band structure that would indicate the softening of these phonons. For 
small c/a ratio the boron chains along the c-directions become unstable against transverse 
vibrations. This is also found to occur for the c/a ratio of 0.675 for lattice parameters 
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FIG. 5: FP-LMTO band structure of bet B for lattice parameter a=4.35 a.u., c/a=0.675 (a), 
a=4.0 a.u., c/a=0.675 (b), and a=4.35 a.u., c/a= 0.65 (c). The zero of energy has been set at the 
Fermi level. 

in excess of a=4.35 a.u. At such large volumes boron presumably enters the nonmetallic 
icosahedral phase. The bands at the Fermi level are mostly of "Px and "Py character. 



B. Phonons and Electron-phonon Coupling 

In FigElthe phonon spectra of bet B for two different lattice parameters and with the c/a 
ratio of 0.675 are shown. The linear response calculations for the phonon properties were 
carried out for 30 points in the IBZ resulting from a 6,6,6 division of the BZ. The solid lines 
in FigEl have been drawn through the calculated frequencies to provide a guide to the eye 
and should not be interpreted as the actual phonon branches. 

With increasing values of the lattice parameter phonons at other symmetry points become 
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FIG. 6: Calculated phonon spectra for bet B for the lattice parameters (a) a=4.25 a.u. and (b) 
a=4.35 a.u. for the c/a ratio of 0.675. The dots represent the calculated frequencies, while the 
solid lines have been drawn through them to guide the eye. For lattice parameters higher than 
4.35 a.u. the structure becomes unstable with the N-point transverse phonons becoming soft first. 

soft as well. Thus our L-R calculations yield real phonon frequencies for a small range of 
volume for the chosen c/a ratio of 0.675. For c/a=0.65 the N-point transverse phonons are 
found to be imaginary for all lattice parameters corresponding to volumes where the bet 
phase should be stable according to the energy-volume curve in Fig^ Of all the phonons 
the ones at the N point are found to be most sensitive to the c/a ratio. 

In Figd the phonon density of states and the Eliashberg spectral functions for three 
different lattice parameters are shown. Unlike in the fee phase the bulk of the electron- 
phonon coupling is via lower frequency phonons, as most of the spectral weight in Fig. 7(b) 
comes from the lower half of the allowed frequency range. As in the fee case, the Eliashberg 
spectral function and the phonon density of states follow each other closely. There is no 
disproportionately large contribution from a particular mode. 

V. SUPERCONDUCTING TRANSITION TEMPERATURES AND SUMMARY 
OF RESULTS 

Superconducting transition temperatures were calculated by solving the isotropic 
Eliashberg equation. We used the calculated Eliashberg spectral function along with the 
following procedure to determine the Coulomb pseudopotential fi*{uc). We start by assum- 
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FIG. 7: Eliashberg spectral function (a) and phonon density of states (b) for bet B, for three 
different lattice parameters, and c/a=0.675. 

ing ^{Ep) = 1.0. This value is consistent with that for MgB2, where the measured ~ 40 
K can be reproduced with the linear response results of phonon frequencies and Eliashberg 
function with /i ~ 0.85. ^*{ujc) was calculated from 

^ ^"""^ 1 + ^^{Ep) HEp/u,) ■ 

Thus from the calculated Fermi energies Ep we obtain fi* for all volumes, with the cut-off 
frequency Uc assumed to be ten times the maximum phonon frequency. A different value of 
^{Ep), e.g. 0.8 or 0.7, would result in a minor change in the calculated and would have 
no effect on the nature of variation of T^. with volume or pressure. 

The results of our calculation are summarized in Table |1J where, for comparison, results 
for MgB2 and two typical low superconductors fee Pb and bcc Nb obtained via the same 
method are also presented. Both for fee and bet boron the Hopfield parameter rj increases 
with decreasing volume per atom. This trend is also supported by the rigid muffin-tin 
approximation results of Papaconstantopoulos and Mehl^l. The Hopfield parameter in these 
high pressure phases of boron are 2-4 times larger than in MgB2 and about 10-15 times larger 
than in typical low metals becauses of the vastly reduced distances between the atoms. 
In fee boron rj is larger than in the bet phase, since the fee phase is the equilibrium phase 
at higher pressures or lower volumes. Phonon frequencies uj in the high pressure bet and fee 
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FIG. 8: Variation of the superconducting transition temperature Tc with volume per atom. The 
vertical dashed line is drawn at volume per atom equal to 21.3 bohr^, estimated to be the dividing 
line between the bet and the fee phases. 



phases of boron are much higher than in MgB2 or a typical metal. As a result the increase 
in the electron-phonon coupling parameter A, which is proportional to r]/ {uj'^) is less drastic 
(than in 77) when compared with MgB2 or a low T^. superconductor like Pb or Nb. The 
area under the Eliashberg spectral function vs. uj curve A is 5-15 times larger for fee or bet 
boron than for a typical electron-phonon superconductor, because the coupling of electrons 
in these high pressure phases takes place via very high frequency phonons. 

The variation of the superconducting transition temperature Tc with volume per atom is 
shown in FiglHl where the vertical dashed line is drawn at volume per atom equal to 21.3 
bohr^, estimated to be the dividing line between the bet and the fee phase. 

For both bet and fee boron the phonon frequencies rise faster with increasing pressure 
than the Hopfield parameter, resulting in a decrease in electron-phonon coupling and Tc. 
In the study of the fee phase by Papaconstantopoulos and Mehlii no calculation of phonon 
spectra as a function of pressure was done, and the authors speculated that if the phonon 
frequencies do not rise as fast as the Hopfield parameter, T^. should rise with increasing 
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pressure. The present calculations reveal two important results: Tc decreases with increasing 
pressure despite rise in the Hopfield parameter and that the calculated Tc's are in the range 
50-100 K, much higher than the measured values. In the fee phase electron-phonon coupling 
rises strongly due to softening of X-point transverse phonons as volume per atom approaches 
the fee— s>bct phase transition. Similarly in the bet phase electron-phonon coupling is very 
high near the bet— s>icosahedral phase change due to softening of N-point transverse phonons. 
Unlike in MgB2 there is no disproportionately high contribution to electron-phonon coupling 
from a particular phonon mode. The Eliashberg spectral function follows faithfully the 
phonon density of states. 

The striking differences between the results of this theoretical investigation and the ex- 
perimental results beg for a suitable explanation. While we cannot rule out the possibility 
that the electron-phonon coupling and therefore Tc is somehow overestimated in the calcu- 
lation, it would be hard to deny the validity of the results altogether. The values of the 
Hopfield parameter and the trend in their variation with pressure are similar to those found 
via a less rigorous rigid muffin-tin approximation calculation of Papaconstantopoulos and 
Mehl^. The phonon calculations correctly predict the fcc-^bct transition and also the onset 
of instability of the bet phase (presumably against the icosahedral phase) as the pressure 
is lowered or the volume per atom is increased. Since the systems studied do not represent 
cases of extreme electron-phonon coupling (see e.g., J. M. An et al.^) via a single mode or 
a few modes, superconductivity in these systems does not seem to be special. 

On the experimental front, the possibility that the measurements were not carried out for 
a single (pure) phase cannot be ruled out, as was pointed out by the authors themselves^. It 
is possible that the systems studied were in mixed phases with co-existing icosahedral and 
bet clusters and the superconductivity of the sample was due to the proximity effect of the 
bet clusters. This would explain why the measured was so much lower than what would be 
expected for the pure bet phase. It is also possible that with increasing pressure the samples 
were just evolving towards pure bet phase, which would explain the observed increase in 
Tc with pressure. The results of the present work strongly suggest that the experimental 
work on high pressure phases of boron should be repeated, with particular emphasis on the 
phase of the sample. It is expected that if the phase is bet or fee (or a mixture of the two), 
superconducting transition temperatures should be much higher than previously reported. 
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TABLE I: Variation of superconducting transition temperature with volume in the bet and fee 
phases of boron. V = volume per atom, to stands for phonon frequency and (• • • )s denote averages, 
Lj = \J (ufi). LOm is the maximum phonon frequency, rj is the Hopfield parameter and A denotes the 
area under the phonon frequency versus Eliashberg spectral function curve a'^F. For comparison, 
results for hexagonal MgB2 and two typical low Tc superconductors, fee Pb and bcc Nb, are also 
presented, denotes the bohr radius. Results for MgB2 are taken from Reference Isl 
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